Breit-Wigner width and participation ratio in finite interacting Fermi systems 
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For many-body Fermi systems we determine the dependence of the Breit-Wigner width V and 
inverse participation ratio £ on interaction strength U > U c and energy excitation 5E > 5E c h when 
a crossover from Poisson to Wigner-Dyson P(s)-statistics takes place. At U > U c the eigenstates 
are composed of a large number of noninteracting states and even for U < U c there is a regime 
where P(s) is close to the Poisson distribution but £ 2> 1. 
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In 1955 Wigner introduced the local density of 
states to study "the properties of the wave functions of 
quantum mechanical systems which are assumed to be 
so complicated that statistical considerations can be ap- 
plied to them". This quantity pw{E) characterizes the 
spreading of eigenstates over the levels of an unperturbed 
system (for example in the absence of interaction between 
particles), and allows to estimate how many of these un- 
perturbed states contribute to the real wave function. 
Generally pw(E) has a Breit-Wigner distribution with 
Lorentzian shape of width T which determines the en- 
ergy spreading over unperturbed states. This concept 
has been shown since then to be very important in a 
wide range of physical problems, from nuclear physics 
and many-electron atoms and molecules to condensed 
matter. 

The study of such complex systems has been success- 
fully performed through the theory of random matrices 
(see for example J|). Very often the physics of such 
systems determines some preferential basis in which the 
Hamiltonian matrix has large diagonal matrix elements, 
while the non-diagonal elements corresponding to tran- 
sitions between the basis states are relatively small. The 
investigation of random matrices of this type has been 
started only recently It has been shown that the 

eigenstates of such superimposed band random matrices 
(SBRM) are spread over the basis states according to 
the Breit-Wigner distribution ; this has been also con- 
firmed analytically through the supersymmetry approach 
[0,1). This spreading determines the number of unper- 
turbed states contributing to a given eigenstate, which 
can be measured through the inverse participation ratio 
(IPR) £ |f-f§- In particular the width T gives an energy 
scale at which the level statistics, for example the num- 
ber variance S2(E), changes behavior from the Wigner- 
Dyson to the Poisson case [||. It has been also shown 
that the Breit-Wigner distribution appears in the case of 
sparse random matrices with preferential basis ||lo| . 

While the properties of the Breit-Wigner distribution 
are well understood in random matrix models, the prob- 
lem of real interacting finite many-body fermionic sys- 
tems was much less investigated. Indeed, in the lat- 



ter case the nature of the two-body interaction should 
be taken into account, since it gives certain restrictions 
on the structure of matrix elements. A very convenient 
model to investigate this kind of problem has been in- 
troduced some time ago in Ojl^). This model consists 
of n fermions distributed over m energy orbitals, cou- 
pled by a random two-body interaction. Recently this 
model attracted a renewed attention since it was under- 
stood that it correctly describes the statistical properties 
of real physical systems such as the Ce atom and the 
Si nucleus p^JI^ |. One of the main advantages of this 
model is that it takes into account the two-body nature 
of the interaction and allows to investigate the depen- 
dence of various quantities on the interaction strength 
U. This property is rather important since the variation 
of the Breit-Wigner width with respect to U and exci- 
tation energy SE counted from the Fermi level has not 
been yet clearly understood. Indeed, due to the two- 
body nature of the interaction, only a small fraction of 
the multiparticle states is coupled by direct transitions. 
As a result, contrary to the common lore the ex- 

ponential growth of the multiparticle density of states p n 
with the number of particles n and the excitation energy 
SE does not imply that an exponentially small interac- 
tion leads to level mixing In a similar way this 
exponential growth of p n does not lead to an exponen- 
tial growth of the width T. This fact has been known in 
nuclear physics for some time |17],0]; however, the precise 
dependence of T on SE has not been determined up to 
now. The dependence of T on U is also not obvious, due 
to the absence of direct coupling between the majority of 
the multiparticle states. Different types of power-law de- 
pendence have been recently proposed JlJ,[ll| but a def- 
inite expression for T is still elusive. A similar situation 
exists for the IPR £ in the basis of noninteracting eigen- 
states which has been studied extensively very recently 
|po| HJ. In this paper, we address these problems and 
determine the dependence of T and £ on the parameters 
above. We show that these two quantities are directly 
related. Surprisingly £ can be arbitrarily large at the 
critical interaction strength U c where the crossover 
in level spacing statistics P(s) between the Poisson and 



1 



Wigner-Dyson distributions takes place. 

To investigate these properties we choose the Two- 
Body Random Interaction Model (TBRIM) described 
above and studied recently in Jl^Jl^]. In this model n 
fermions are located on m orbitals with one-particle ener- 
gies e' m randomly chosen in the interval [0, m] so that the 
average one-particle level spacing is A = 1. The multi- 
particle states are distributed from the ground state E g « 
n 2 A/2 to the maximal energy E t w mnA — E g . These 
states are coupled by two-body random matrix elements, 
varying in the interval [-U, U]. Due to the two-body na- 
ture of the interaction, a given multi-particle state is only 
coupled to K = l+n(m— n)+n(n— l)(m— n)(m— n— 1)/4 
other states in an energy interval B = 2m — 4. This 
K is much smaller than the total number of states 
N = m\/n\(rn — n)l. The density of directly coupled 
states p c = K/B w mn 2 /8 is therefore much smaller 
than the total density p n rj N/(E t — E g ). According 
to the results in a crossover for P(s) from Pois- 

son to Wigner-Dyson statistics takes place at a critical 
interaction strength U c — C/p c with C ~ 0.58. A simi- 
lar border was also discussed in p3| . The precise value 
of U c [[l6| was determined by the condition that rj = 
J So (P(s) - P WD (s))ds/ C(P P (s) - P W D(s))ds = 0.3. 
Here Pp(s) and Pwd(s) are the Poisson and the Wigner- 
Dyson distributions respectively and s = 0.4729... is 
their intersection point. Physically this crossover hap- 
pens when the coupling matrix elements become compa- 
rable to the energy spacings between directly coupled 
states |l8| , |l^ , |l(| . A similar condition determines the 
metal-insulator transition in the Anderson model, where 
also the level statistics P(s) changes from the Poisson 
distribution to the Wigner-Dyson one | p5| . However, the 
TBRIM case differs from the Anderson model where at 
large system size r] can take only three values r/ — 1 (lo- 
calized), 77 = 0.215 (critical), rj = (delocalized), while 
in the TBRIM 77 varies smoothly near U c |l6| . Physically 
this difference comes from the fact that in the Anderson 
model the number of coupled neighbors is much smaller 
than the linear system size, while in the TBRIM this 
number is of the order of the number of states m in one 
of n directions associated to each particle. 

While the value of U c has been determined |ll|, the 
properties of eigenstates as a function of the interaction 
remained unclear. To understand these properties in the 
TBRIM we studied the local density of states pw(E) in 
the basis of noninteracting multiparticle states. The data 
were obtained for the states near the middle of the spec- 
trum (±25% from the center). The total statistics for 
pw was kept around 10 6 . We checked that pw{E) has a 
Breit-Wigner shape and analyzed the dependence of its 
width r on U. The numerical data for TBRIM clearly 
demonstrate the relation r oc U 2 which continues up to 
large U values where a saturation takes place (Fig. 1). 
To check this dependence for larger system sizes, we in- 
vestigated a slightly different model, obtained from the 



TBRIM by restricting ourselves to states in an energy 
layer of width A near total energy E = mA. Such an 
approximation is physically reasonable provided r <C A. 
Indeed, in this case the transitions to states outside the 
layer do not influence the properties of eigenstates. We 
choose the layer to be defined by Y^i=i m 'i = m - The 
transition matrix elements between these states were 
taken from the TBRIM, and the diagonal elements com- 
ing from one-particle energies e' m were randomly chosen 
in [(m - 1/2) A, (m + 1/2) A]. The layer model (LM) de- 
fined in this way retains the main physical properties of 
the TBRIM but allows to study systems with much larger 
number of orbitals m. For n — 3 the system size of LM 
is N w to 2 / 12 and for n = 4, N w m 3 /200. This allowed 
to span m-values up to m — 130 (n = 3, N w 3.6 10 5 ) 
and m = 60 (n = 4,A«4.9 10 6 ), which are much larger 
than the values reached in |l^Jl^,|2^] . The multi-particle 
density in the LM is p n — N / A while p c was determined 
numerically. The data for LM (Fig. 1), similarly to the 
TBRIM case, also demonstrate the dependence r cx U 2 
and show in addition that r oc p c . 
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FIG. 1. Dependence of the rescaled Breit-Wigner width 
T/(p c A 2 ) on 17/ A: TBRIM data for n = 3, m = 17 (X); LM 
data for n = 3,m = 130 (o) and n = 4,m = 60 (diamonds). 
The full line shows the theoretical estimate (1). Insert gives 
an example of pw{E) for LM (*) with the Breit-Wigner fit 
(r = 0.12) for n = 3,m = 130, U = 0.022 when (1) gives 
T = 0.125. 

According to the data of Fig. 1 the width T is given 
by the Fermi golden rule: 

2tt, 



T = 2tt < U 2 > p c 



-U 2 Pc . 



(1) 



where < ... > means the averaging. We attribute the 
small difference between the LM and TBRIM cases to the 
fact that in the latter the density p c slightly depends on 
the energy counted from the Fermi level, while we used 
its average value. For the LM this variation is smaller 
and therefore the agreement is better. 
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The expression (|l|) for T doesn't depend on the multi- 
particle density of states p n . However, for U > U c the 
levels are mixed on a scale A„ = l/p„ *C 1//0 C ■ There- 
fore we expect that an eigenstate is spread over all un- 
perturbed states in the energy interval T; for U > U c this 
leads to the IPR: 

£«I>„ « 2U 2 PcPn . (2) 

The numerical factor was taken in analogy with SBRM 
case where £ w Tp [0J|]- To check this theoretical es- 
timate we computed £ for both TBRIM and LM. The 
numerical data displayed in Fig. 2 show clearly the U 2 
dependence for sufficiently strong U. At very large U the 
growth of £ is replaced by a saturation due to the finite 
size of the system. The data shown in Fig. 3 demonstrate 
that £ oc p c Pn, in agreement with (2). 
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FIG. 2. Dependence of the IPR £ on U/A: TBRIM data 
for n — 3,m = 17 (X); LM data for n = 3, m = 130 (o), 
n — 3, m = 90 (triangles) and n = 4, m = 60 (diamonds). 
Straight lines show dependence £ oc U . 

Without any fitting parameters these numerical results 
definitely confirm the estimate (2) for U > U c . It is 
interesting to check if it remains valid close to the critical 
value U c . If so, then the IPR at U c for n ^ 1 will contain 
exponentially many states £ c = = U c ) ~ p n /p c 3> 1. 
We studied this behavior in both TBRIM and LM. In the 
latter case, we checked that U c defined by the condition 
i](U c ) = 0.3 also follows the relation U c — C/p c , with 
C = 0.62 being very close to the TBRIM value |j| (see 
insert in Fig. 3). This fact once more confirms that 
indeed the LM retains the physical properties of TBRIM. 
The data of Fig. 4 for £ c confirm that £ c ~ p n /p c in 
both TBRIM and LM. In LM the proportionality factor 
C is about 3 times smaller (C ~ pdc/Pn ~ 0.25 in Fig. 
4) than its value given by (2) at U = U c (C w 0.8). 
This indicates a change of eigenstate properties near U c . 
The difference of C values for LM and TBRIM should be 
attributed to a stronger variation of the densities p n ,Pc 



with energy in the TBRIM. This variation was not taken 
into account in the expressions for p n ,p c in the TBRIM 
where we used their averaged values. 
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FIG. 3. Dependence of the rescaled IPR £,/U 2 on p c p n : LM 
data for n — 3 and 40 < m < 130 (o); n — 4 and 30 < m < 60 
(diamonds). The straight line gives theory (2). Insert shows 
Uc/A vs. p c A in log-log scale for the same parameters; the 
straight line is the fit U c = 0.62/p c . 

The data of Fig. 4 definitely show that at U = U c 
the IPR grows proportionally to the multi-particle den- 
sity p n and, therefore, it is exponentially large for large 
number of particles n. This fact leads to the apparently 
surprising conclusion that for U < U c the eigenstates are 
composed of a huge number of noninteracting eigenstates 
but the level statistics P(s) is still close to the Poisson 
distribution. A similar situation is known to exist for 
quantum systems whose classical dynamics corresponds 
to the Kolmogorov - Arnold - Moser (KAM) regime. In 
this case, the coupling between different modes strongly 
deforms the unperturbed tori, but doesn't destroys the 
integrals of motion and the corresponding quantum num- 
bers. Generally such deformation gives a spreading over 
many unperturbed eigenstates, without real mixing of 
energy levels [^6| . The mixing and Wigner-Dyson statis- 
tics for -P(s) appear only after the transition to chaos, 
which in our case corresponds to the situation when the 
physical frequency l/p c becomes comparable to the in- 
teraction induced interstate transition rate T. A similar 
phenomenon takes place near the Anderson transition 
|p5| . However, in this case £ becomes infinite at/above 
the transition. 

We have so far discussed the case of highly excited 
states far from the Fermi level, where p n and p c are not 
very sensitive to energy variation. This is not so near 
the Fermi energy ep ~ nA, where the dependence on 
excitation energy SE = E — E g should be taken into 
account. At temperature T only Sn ~ T/A particles 
effectively interact near the Fermi level so that SE ~ 
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TSn ~ T 2 /A. Since near tp the density p c ~ piiSn) 
(Sn) 3 /A 0, we obtain 

C/ 2 /<5£\ 3/2 £/ 2 /T\ 2 r 



A~ 



1/2 



exp 2 



tt 2 (5£; 

6A 



1/2N 



(3) 



Here we used the known dependence of p„ on SE from 
|3 and assumed that SE > SE ch « A(A/£f) 2/3 (T > 
T c h « A(A/C/) 1 / 3 ) so that the system is thcrmalized due 
to internal interaction (U > U c ) [^6| . The last expression 
for r has a simple meaning. Indeed, T is the total spread 
width for Sn effectively interacting particles. Therefore, 
the partial width ~ T /Sn is the usual quasi-particle 
decay rate which in agreement with the theory of Landau 
Fermi liquid is proportional to T 2 . At the quantum chaos 
border SE = SE^, when the crossover to the Wigner- 
Dyson statistics takes place , the IPR becomes expo- 
nentially large £ c ~ (T ch /A) exp(2.6T c/l /A). 
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FIG. 4. Dependence of the IPR £ c at U c on p„/p c : LM data 
for n = 2 and m = 800 (*); n = 3 and 40 < m < 130 (o); 
n = 4 and 30 < m < 60 (diamonds). The straight line gives 
the fit £ c = 0.25p„/p c + 2.7. Insert shows the same plot for 
the TBRIM for n = 2, m = 30; 3 < n < 6 and 10 < m < 21 
(*). The straight line is the fit £ c = 1.63p„/p c + 1.91. 

In conclusion, our results allowed to understand the 
eigenstate properties in finite Fermi systems with in- 
terparticle interaction U > U c ~ l/p c and excitation 
energy SE > 5E c h- Qualitatively, these properties are 
similar and opposite to the recent expectations p2| and 
]l~3| p^| , p^| pl| ] correspondingly. Further investigations are 
required for the regime with U < U c where there are indi- 
cations on the appearance of another dependence of the 
IPR on system parameters |^3| . Direct numerical checks 
of the relations (3) are also desirable. 

We thank Ph.Jacquod for useful discussions. 
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